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Abstract 

Several new effects have been investigated in recent analyses of supersymmetric 
dark matter. These include the effects of the uncertainties of wimp velocity dis- 
tributions, of the uncertainties of quark densities, of large CP violating phases, of 
nonuniversalities of the soft SUSY breaking parameters at the unification scale and 
of coannihilation on supersymmetric dark matter. We review here some of these with 
emphasis on the effects of nonuniversalities of the gaugino masses at the unification 
scale on the neutralino-proton cross-section from scalar interactions. The review 
encompasses several models where gaugino mass nonuniversalities occur including 
SUGRA models and D brane models. One finds that gaugino mass nonuniversalities 
can increase the scalar cross-sections by as much as a factor of 10 and also signifi- 
cantly extend the allowed range of the neutralino mass consistent with constraints 
up to about 500 GeV. These results have important implications for the search for 
supersymmetric dark matter. 



1 INTRODUCTION 



Over the recent past there has been considerable experimental activity in the di- 
rect detection of dark matter E> @) and further progress is expected in the ongoing 
experiments 

I, I 1) 

and new experiments that may come online in the future o) . 
At the same time there have been several theoretical developments which have shed 
light on the ambiguities and possible corrections that might be associated with the 
predictions on supersymmetric dark matter. These consist of the effects on the dark 
matter analyses of wimp velocity 0' 0' 0) and of the rotation of the galaxy ID , the 
effects of the uncertainties of quark densities and the uncertainties of 

the SUSY parameters 0), effects of large CP violating phases 0> 0), effects of 
scalar nonuniversalities 0) , effects of nonuniversalities of gaugino masses LLJJ/ and 



effects of coannihilation 0) . In this paper we will discuss some of these briefly but 
mainly focus on the effects of nonuniversalities of the gaugino masses on dark matter. 
In the Minimal Supersymmetric Standard Model (MSSM) there are 32 supersym- 
metric particles and with R parity conservation the lowest mass supersymmetric 
particle (LSP) is absolutely stable. In many unified models, such as in the SUGRA 
models 0), one finds that the lightest neutralino is the LSP over most of the pa- 
rameter space of the model. Thus the lightest neutralino is a candidate for cold 
dark matter. The quantity that constrains supersymmetric models is Q x h 2 where 
Q x = p x j p c where p x is the density of relic neutralinos at the current temperatures, 
and p c = 3Hq/8tcGn is the critical matter density, and h is the Hubble parameter 
Hq in units of 100 km/sMpc. The most recent measurements of h from the Hubble 
Space Telescope give 0) 

h = 0.71 ±0.03 ±0.07 (1) 
Similarly the most recent analyses of Vl m give 0) 

tt m = 0.3 ± 0.08 (2) 

If we assume that the component of f2g m is — 0.05 which appears reasonable, 
then this leads to the result Q x h 2 = 0.126 ± 0.043. Perhaps a more cautious choice 
of the range would be a ~ 2a range which gives 

0.02 < tt x h 2 < 0.3 (3) 



The quantity of interest theoretically is 



A/2 



n v ^ 2 .48 Xl o-"(!)W^ < 4 > 



Here Tf is the freeze-out temperature, Xf = kTf/m^ where k is the Boltzman 
constant, Nj is the number of degrees of freedom at the time of the freeze-out, (p^) 3 
is the reheating factor, and J(xf) is given by 

J (x f ) = P dx (av) {x)GeV~ 2 (5) 
Jo 

where < av > is the thermal average with a the neutralino annihilation cross-section 
and v the neutralino relative velocity. 

2 Detection of Milky Way wimps 

Both direct and indirect methods are desirable and complementary for the detec- 
tion of Milky Way wimps. We shall focus here on the direct detection. In this 
case the fundamental detector is the quark and the relevant interactions are the su- 
pergravity neutralino-quark-squark interactions. The scattering of neutralinos from 
quarks contains squark poles in the s channel and the Z boson and the Higgs bo- 
son (h, H°, A ) poles in the t channel. Since the wimp scattering from quarks is 
occuring at rather low energies one may, to a good approximation, integrate on the 
intermediate squark, Z and Higgs poles to obtain a low energy effective Lagrangian 
which gives a four-Fermi interaction of the following form 

£eff = XIh^XQ^^Pl + BP R )q + Cxx™qM + #X75X""W5<7 (6) 

The contribution of D is generally small and thus the scattering is effectively gov- 
erned by the terms A,B and C. Analysis of dark matter is affected by several factors. 
We discuss these briefly below. 

2.1 Uncertainties in wimp density and velocity 

Two of the quantities that control the detection of dark matter are the wimp mass 
density and the wimp velocity. Estimates of Milky Way wimp density lie in the 
range p x = (0.2 — 0.7)GeVcm~ 3 and the event rates in the direct detection 
depend directly on this density. A second important factor regarding wimps that 
enters in the dark matter analyses is the wimp velocity. One typically assumes a 
Maxwellian velocity distribution for the wimps and the current estimates for the rms 
wimp velocity give v = 270 km/s with, however, a significant uncertainty. Estimates 
for the uncertainty lie in the range of ±24 km/s to ±70 km/s 0). A reasonable 
estimate then is that the rms wimp velocity lies in the range 0) 



v = 270 ± 50 km/s 



(7) 



Analyses including the wimp velocity variations show that the detection rates can 
have a significant variation, i.e., a factor of 2-3 on either side of the central values 0). 



2.2 Effects of uncertainties of quark densities 

The scattering of neutralinos from quarks are dominated by the scalar interaction 
which is controlled by the term C in Eq.(6). The dominant part of the scattering 
thus arises from the scalar part of the x~ P cross-section which is given by 

Au 2 2 
a XP (scalar) = -A( £ ff Q + -(1 - £ ff) £ C a f (8) 

i=u,d,s i=u,d,s a=c,b,t 

Here /i r is the reduced mass in the x~ P system and ff (i=u,d,s quarks) are quark 
densities inside the proton defined by 

m pfi =< p\m qi qiqi\p > (9) 

There are significant uncertainties in the determination of ff. To see the range of 
these uncertainties it is useful to parametrize the quark densities so that 0) 
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f!= ==—(1-0;)^ (10) 



where we have defined 



o"o < p\uu + dd — 2ss\p > <p\uu — dd\p> 

X = = — F = — 

O'-kN < p\uu + dd\p > < p\uu + dd\p > 

ct^n =< p\2~ 1 (m u + m d )(uu + dd\p > (11) 
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The current range of determinations of u^n give 

a nN = 48 ± 9 MeV, x = 0.74 ± 0.25, ^ = 0.132 ± 0.035 (12) 

With the above range of errors one finds that ff lie in the range 

ff = 0.021 ±0.004, ff = 0.029 ±0.006, ff = 0.21 ± 0.12 (13) 

Of these the errors in ff generates the largest variations. A detailed analysis shows 
that the scalar cross-section can vary by a factor of 5 in either direction due to errors 
in the quark densities O) . 



2.3 CP violation effects on dark matter 



The soft SUSY breaking parameters that arise in supersymmetric theories after 
spontaneous breaking are in general complex with phases 0(1) and can lead to large 
electric dipole moment of the electron and of the neutron in conflict with current 
experiment. Recently a cancellation mechanism was proposed as a possible solution 
to this problem H3' 2D . With the cancellation mechanism the total EDM of the 
electron and of the neutron can be in conformity with data even with phases 0(1) 
and sparticle masses which are relatively light. The presence of large CP phases can 
affect dark matter and other low energy phenomena. Effects of CP violating phases 
on dark matter have been investigated for some time 2D and more recently such 
analyses have been extended to determine the effects of large CP phases under the 
cancellation mechanism 0) . One finds that the EDM constraints play a crucial role 
in these analyses. Thus in the absence of CP violating phases one finds that the 
X — p cross-section can change by orders of magnitude when plotted as a function 
of 6^ (the phase of //). These effects are, however, significantly reduced when the 
constraints arising from the current experimental limits on the electron and on the 
neutron EDMs are imposed. After the imposition of the constraints the effects of 
CP violating phases are still quite significant in that the x~P cross-section can vary 
by a factor of ~ 2. Thus precision predictions of the x~P cross-section should take 
account of the CP phases if indeed such phases do exist in a given model. Indeed 
many string and D brane models do indeed possess such phases and thus inclusion of 
such phases is imperative in making predictions for direct detection in such models. 



2.4 Effects of coannihilation 

The effects of coannihilation may become important when the next to the lowest 
supersymmetric particle (NLSP) has a mass which lies close to the LSP mass 0). 
The size of the effects is exponentially damped by the factor e~ AlX where Aj = 
(mi/m x — 1), x = m x /kT and where m x is the LSP mass. Because of this damping 
the coannihilation effects are typically important only for regions of the parameter 
space where the constraint Aj < 0.1 is satisfied. Some of the possible candidates for 
NLSP are the light stau fi, Br, the next to the lightest neutralino x% an d the light 
chargino xt ■ An interesting result one finds is that in mSUGRA the upper limit 
on the neutralino mass consistent with the current experimental constraints on the 
relic density is extended from 200 GeV to 600 GeV 0) when the effects of x ~ t 
coannihilation are included. In Sees. (2.6) and (2.7) we will show that the allowed 
range of the neutralino mass can also be extended by inclusion of nonuniversalities 



in the gaugino masses. 



2.5 Nonuniversality of scalar masses 

The minimal SUGRA model is based on the universality at the GUT scale. This 
includes the universality of the scalar masses, of the gaugino masses and of the tri- 
linear couplings at the GUT scale. In supergravity unified models the universality 
of the soft SUSY breaking parameters arises from the assumption of a flat Kahler 
potential. However, the nature of physics at the GUT/Planck scale is not fully un- 
derstood and a more general analysis of the soft SUSY breaking sector requires that 
one work with a curved Kahler potential El' HI) . Such an analysis in general leads 
to nonuniversalities in the scalar sector of the theory. However, the nonuniversali- 
ties in the scalar sector cannot be completely arbitrary as there are very stringent 
constraints on the system from the limits on the flavor changing neutral currents 
(FCNC). A satisfaction of the constraints requires essentially a degeneracy in the 
scalar masses in the first two generations at the GUT scale. However, the constraints 
on the scalar masses in the Higgs sector and on masses in the third generation are 
far less severe and one could introduce significant amounts of nonuniversalities in 
these sectors without violating the FCNC constraints. It is found convenient to 
parametrize the nonuniversalities in the Higgs sector by 81,82 so that at the GUT 
scale (M G ) one has m 2 H = wig(l + 8\), m^ 2 = m§(l + <5 2 ). Similarly one may 
parametrize the nonuniversalities in the third generation squark sector by £3, £4 so 
that at the scale Mc one has m 2 ^ = 77^(1 + £3), mf r = m\ (1 + 84). These nonuni- 
versalities have a significant effect on the low energy physics. One of the main effects 
that occurs is through the effect on \i which is determined via the constraint of the 
radiative breaking of the electro-weak symmetry and is modified in the presence of 
the nonuniversalities in the Higgs sector and in the third generation sector. To one 
loop order it is given by 

f = nl + ^(81 - 8 2 t 2 - \±{8 2 + 5 3 + 5 4 )t 2 ) + Afi 2 (14) 

Here /x is the value of fi in the absence of nonuniversalities, D depends on the 
top Yukawa coupling and defines the position of the Landau pole, t = tan/3, and 
Afj 2 is the loop correction. We note that the entire effect of nonuniversalities is now 
explicity exhibited. One finds that the universalities can significantly affect the event 
rates. The effect on the event rates occurs specifically because of the effect on fi. 
Thus one finds that for certain regions of the parameter space the nonuniversalities 
in the Higgs and in the third generation sector make a negative contribution to fi 2 



which leads to larger higgsino components for the neutralino. Since in the direct 
detection the scattering is dominated by the scalar x~P cross-section which in turn 
depends on the product of the gaugino and the higgsino components one finds that 
a smaller /i leads to larger event rates in the direct detection. A detailed analysis 
of the effects of nonuniversalities of the scalar masses has been given in Refs. . 
We will discuss further this phenomena in the context of the nonuniversalities in the 
gaugino sector in the next section. 



2.6 Gaugino nonuniversalities and dark matter in GUT models 

Nonuniversality of gaugino masses arises in grand unified models via corrections to 
the gauge kinetic energy functions cM) . Thus in grand unified models a non-trivial 
gauge kinetic energy function leads to a gaugino mass matrix which has the form HI) 

map = \e G/2 G a {G^) b a {dr ai /dz* b )f^ (15) 

As an example if we consider the GUT group to be SU(5) then the gauge kinetic 
energy function f a p transforms as follows 



f24 x 24 



symm 



1 + 24 + 75 + 200 



(16) 



where ( 24 x 24' 



symm 



stands for the symmetric product. The term that transforms 



like the singlet of SU(5) in the gauge kinetic energy function leads to universality 
of the gaugino masses, while the 24 plet, the 75 plet and the 120 plet will generate 
corrections to universality. In general one could have an admixture of the various 
representations and this will lead to gaugino masses of the form 

!7w(0)=mi(l + 5>n*) (17) 



where n\ depend on r and for the representations 1, 24, 75, 210 they are given in 
Table 1 . The nonuniversality of the gaugino masses also leads to corrections of 
the gauge coupling constants at the GUT scale and in general one has q%{Mg) = 
9g(1 + Sr c r n D- We note, however, that the coefficients c' r that enter in ^ are 
different than those that enter in mj. This is so because the corrections to g± involve 
only the gauge kinetic energy function while the corrections to mi involve the gauge 
kinetic energy function as well as the nature of GUT physics. 

Table: nonuniversalities at Mx- 



SU(5) rep 


n\ 


n r 2 


T 

n 3 


1 


1 


1 


1 


24 


-1 


-3 


2 


75 


-5 


3 


1 


200 


10 


2 


1 



The gaugino sector nonuniversalities affect \i. To exhibit this effect we can expand \i 
determined via the constraint of the radiative breaking of the electro-weak symmetry 
in terms of the parameter c r . One finds the following expansion 

A^o + Zf^V + O^) (18) 



and for c 2 4 < 0, C75 < 0, C200 > one has 
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^li >0 ,^>0,^<0 (19) 
dc 2A oc 75 dc 200 

Thus in these cases the nonuniversalities lead to a smaller value of Now as 
already mentioned in the previous section the Higgsino components become more 
dominant as /i becomes smaller. We can exhibit this analytically for the case when 



H is small but we are still in the scaling region where /i 2 /M| >> 1. In this case 



it is possible to analytically investigate the size of the gaugino- Higgsino components 
X n0 of the LSP defined by 

X = X 10 B + X 20 W 3 + X 30 #i + X m E 2 (20) 

where B is the Bino, W3 is the Wino, and Hi, and H 2 are the two Higgsinos. In 
this case one finds that the gaugino components of the LSP are given by Xu ~ 
1 — (^-r)sin 2 6 w , and X 12 — , M 2 Z sin26wsin(3 while the higgsino components are 

given by il) X13 ~ ——sin29w sin/3, Xu ~ — sin29w sinf3 . From the above one 
finds that the Higgsino components have a dependence on the inverse power of /i and 
thus a smaller fi will lead to a larger scalar u x - p cross-section. The literature on the 
analyses of dark matter relic density 0) and direct detection in MSSM 0) and 
in SUGRA models 0) is quite extensive © . We discuss here the quantitative effects 
of the gaugino mass nonuniversality on dark matter. Some features of the effects 



of gaugino mass nonuniversalities have already been discussed in the literature 0) 
and we review here the more recent developments O) . The techniques used in the 
analysis are as discussed in Ref. 0) and in the analysis we impose the 6^5 + 7 
constraint HI) . In Fig.|I] we plot the scalar x ~ P cross-section as a function of the 
neutralino mass for the case of GUT scale nonuniversalities with values of C24 in the 
range -0.1 to 0.08. One finds that the scalar cross-section is enhanced for negative 
values of c 2 4 just as one would expect from the general discussion above because 
it is for the case of c 24 negative that fi becomes small. One finds that in general 
the scalar cross-section increases systematically as \c 2 ^\ increases for negative values 
of c 2 4 and an enhancement of the scalar cross-section by as much as a factor of 10 
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Figure 1: a xp (scalar) vs m x when m = 51 GeV, tan/3 = 10, A t /m = —7 and c 2 4 
takes on various values (Taken from Ref. H3) ). 

can be gotten relative to the universal case of C24 = 0. One also finds an enhance- 
ment of the allowed range of the neutralino mass consistent with the constraints. 
In Fig.^| we plot the maximum and the minimum of the scalar cross-section as a 
function of the neutralino mass for the case of GUT scale nonuniversalities where 
the nonuniversalities arise from the 200 plet representation with C200 = 0.1 when 
the other parameters are varied over their assumed naturalness range. The current 
experimental limits from DAMA E) and from CDMS ED are also plotted. Further, 
the currents limits would certainly be significantly improved in other dark matter 
detectors in the future H' & H) and in Fig.^| we also plot the expected limits from 
future CDMS, and from GENIUS ®. One finds that the current experiment does 
constrain the theory in a small region of the parameter space. Further, the expected 
sensitivity in future experiment, i.e., in CDMS and in GENIUS will explore a ma- 
jor part of the parameter space of this model. We also note that the inclusion of 
nonuniversality significantly increases the allowed range of the neutralino parameter 
space. 

2.7 Dark Matter on D Branes 

Nonuniversality of gaugino masses is rather generic in string theory. However, the 
specific nature of the nonuniversality will depend on the details of the compactifi- 
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Figure 2: The maximum and the minimum curves of a xp (scalar) vs m x when 
C200 — 0.1 and SUSY parameters are varied over the natualness domain (Taken 

from Ref. HI) ). 



cation. We discuss here the effects of gaugino mass nonuniversality on dark matter 
in the context of D brane models. The possibility of nonuniversal gaugino phases in 
brane models arises from the choice of embedding of the different gauge groups of the 
Standard Model on different branes. One may consider, from example, models that 
arise from Type IIB string compactified on a six-torus T 2 x T 2 x T 2 which contains 
9 branes, 7j and 5j (i=l,2,3) branes and 3 branes. Not all the branes can be present 
simultaneuosly due to the constraint of N=l supersymmetry which requires that 
one has either 9 branes and 5j (i=l,2,3) branes or 7j (i=l,2,3) branes and 3 branes. 



In the following we will make the choice of embedding on 9 branes and 5 branes 
One of the major problems in developing a sensible string phenomenology is that 
the mechanism of supersymmetry breaking in string theory is still lacking. However, 
some progress can be made by use of an efficient parametrization of supersymmetry 
breaking. Here we use the parametrization where the breaking of supersymmetry 
arises from the breaking generated by the dilaton and the moduli VEV's of the fol- 
lowing form ED F s =yj3m 3/2 (S + S*)sin9e^ s , F i =V^m 3/2 (T i + T^cosOie'^ . 
We consider now a specific 9-5 brane model. Here one embeds the SU(3)c x U(\)y 
gauge group on 9 branes and SU(2)l gauge group on a 5i brane. The alternative 
possibility of embedding the Standard Model gauge group on five branes is discussed 
in the last two papers of Ref. EH) . For the 9 — 5i brane model the soft SUSY breaking 
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Figure 3: The maximum and the minimum curves of a xp (scalar) vs m x for the 9-5 
D brane model when C200 — 0.1 when m 3 / 2 ranges up to 2 TeV, tan/3 ranges up to 

25, and 9 lies in the range 0.1-1.6. (Taken from Ref. COP ). 
sector of the theory is given by © 0) 



ml 



mi — 1^3 

m 3/2(! 



V3m^/2sin8e ns 
3cos 2 9®\); 



.2 



m 



3/2 



rri2 = V3iri3/2COs8e %l1 
l - (3/2)cos 2 9(l - 



(21) 



Here 9(Qi) is the Goldstino direction in the dilaton S (moduli Tj) VEV space. We 
discuss now dark matter on D branes. In Fig|| we give a plot of the scalar \ ~ P cross 
section as a function of the neutralino mass for the 9 — 5i D brane model. One of 
the interesting feature of the D brane model is that the scalar masses are in general 
not universal. However for Oi = 1/a/3 one has m 9 = and the scalar masses are 
universal although the gaugino masses are still nonuniversal. Since we are mostly 
interested here in investigating the effects of nonuniversalities of the gaugino masses 
in this analysis, we impose universality of the scalar masses and set Q% = l/y/3. In 
Fig.|3| we give a plot of the minimum and the maximum of the scalar \ — p cross 
section under this constraint. One finds that under the assumed constraints the 
allowed domain of the parameter space has the general features which are similar to 
the GUT scale nonuniversalities. One common feature is that the allowed domain of 
the parameter space is extended close to 500 GeV. One may note that if in addition 
to the constraint 61 = one also sets 9 = 7r/6 one finds also universality of the 

gaugino masses. This situation is exhibited by the vertical dark line in the enclosed 



region on the left hand side in Fig.|. 



3 Conclusion 



In this paper we have given a brief review of the recent theoretical developments in 
the analyses of supersymmeteric dark matter. Our emphasis has been in exploring 
the effects of uncertainties of the input data and the effects of nonuniversalities of 
the gaugino masses on dark matter analyses. It is found that the uncertainties of 
the wimp velocities can change detection rates by up to factors of 2-3 while the 
uncertainties in quark masses and densities can change the x ~ P cross-section by 
up to factors of 5 in either direction. The effects of gaugino mass nonuniversalities 
on dark matter analyses is found to be quite dramatic. It is seen that gaugino mass 
nonuniversalities can increase the Higgsino components of the LSP and significantly 
increase the x~P cross-section from scalar interactions and also increase the allowed 
range of the LSP consistent with relic density constraints. Thus an increase in the 
scalar x ~ P cross-sections by up to a factor of 10 can occur while the allowed 
range of the neutralino masses can move up to 500 GeV consistent with the relic 
density constraints. Data from current dark matter experiments is beginning to put 
constraints on models with nonuniversalities. These constraints will become more 
severe as the sensitivity of dark matter experiments increase in the future. 
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